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We study center vortex free energies and 't Hooft's electric fluxes on the lattice in 2+1 dimen- 
sions, where SU{2) for example, is in the universality class of the 2d Ising model. This places 
a wealth of exact results at our fingertips. In particular, spacelike center vortices in SU (2) near 
criticaUty correspond to spin interfaces in the 2d Ising model, whose universal scaling functions 
are known exactly. We exploit this to locate the deconfinement transition with unprecedented 
precision and subsequently for a finite size scaling analysis, where the self-duality of the 2d spin 
model is reflected in a duality between the spacelike vortices and confining electric fluxes. The 
corresponding relation between the string tension and its dual in the high temperature phase is 
arguably the simplest example of a universal amplitude ratio. Around the transition, both can be 
efficiently extracted from the exact results with a global one-parameter fit which allows straight- 
forward continuum extrapolation. 
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1. Introduction 

The finite temperature deconfinement transition in pure SU {N) gauge tlieories in J + 1 dimen- 
sions is very well understood in terms of the spontaneous breakdown of their global center 
symmetry. This symmetry is faithfully represented by the fundamental Polakov loops which live 
on d dimensions and develop a non-zero expectation value in the deconfined broken phase much 
Uke the spontaneous magnetization in a J dimensional g-state Potts model [1] with q = N. 

In order to be able to apply the powerful tools of universality and scaling near a critical point, 
here we are particularly interested in cases where this transition is of 2"'^ order. For pure QCD with 
N = 3 colors in 3 + 1 dimensions this is of course not the case. The transition is first order, but only 
just. If we either reduce the number of colors to A'^ = 2 or the dimensions to 2 + 1 , we obtain Yang- 
Mills theories with a second order deconfinement transition within the universality class of <7-state 
Potts models with q = 2 (Ising) in J = 3 dimensions or with q = 2 and 3 in J = 2. The q = 4 Potts 
model in 2 dimensions is interesting because it is known [2] to have a 2"*^ order ti'ansition and to 
fall precisely on the separatrix qc{d) with respect to the order of the transition in the (^,fif)-plane, 
where only 1**' order transitions occur for q > qc{d), i.e., qc{2) = 4. The corresponding SU{A) 
gauge theory in 2-1-1 dimensions has recently been studied for example in [3, 4, 5]. The conclusion 
in [4] is that the transition is weakly P' order, unlike the q = A Potts case. It nevertheless seems that 
there is a rather wide range of intermediate length scales where at least approximate Potts scaling 
can be observed. One might then like to understand, for example, why among the wider class of 
Ashkin-Teller models it is the standai^d q = A Potts scaling that is relevant here, and whether this 
can be derived from an effective Polyakov-loop model. 

We will present more results specific for SU{3) and SU{4-) in 2-h1 dimensions in [6]. One 
general property of the 2d Potts models is that they are self-dual for all q. We will demonstrate 
below, how this self-duality is reflected in a duality between the spacelike center vortices and the 
confining electric fluxes of the gauge theory. Herein we will mainly report on numerical results for 
SU (2) in 2-1-1 dimensions for which many exact results from the 2d Ising model are available. 

2. Center Vortices and Interfaces 

For pure SU (N) gauge theories in a finite Euclidean box, 't Hooft's twisted boundary condi- 
tions fix the total Zat- valued center flux through each plane of the box, i.e., the total number modulo 
A'^ of center vortices piercing that plane [7]. At finite temperature T, the differences in free energy 
between these twisted ensembles and the periodic one define the corresponding center-vortex free 
energies. If they tend to zero in the thermodynamic limit, then these vortices condense which leads 
to an area law for those Wilson loops that feel their disordering phases in the corresponding plane 
of the twist. In an L'' X 1/r box, we have to distinguish between temporal and magnetic twist. 
The latter is defined in purely spatial planes and corresponds to the ZAf-valued magnetic flux m of 
static center monopoles. Their free energy tends to zero with L — )• oo as exp{—as{T)L^} at all T 
[8], corresponding to the area law for spatial Wilson loops with spatial string tension Gs{T). 

Relevant to the deconfinement transition, however, are only the temporal twists. These ai^e 
characterized by cf-dimensional vectors k of integers mod A'^, i.e., with components kj G rep- 
resenting the twist in the temporal plane of orientation (0,/), with total center flux exp{27iiki/N) 
through that plane. See Fig. 1 for an illustration of such a vortex ensemble in SU (2). We denote the 
partition functions of the various ensembles with temporally twisted boundary conditions by Z^ik). 
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Figure 1: Spatial center vortices share their universal behavior with interfaces in a spin model. 



The corresponding center- vortex free energies per T are then given by Fi^(k) = — ln(Zi(^)/Z^(0)), 
where Zk{0) stands for the periodic ensemble. In the vicinity of a second order deconfinement 
transition these vortex free energies show the universal behavior of interfaces in the respective 
d dimensional Potts model. These interface free energies, Fj{c) = — ln(Z^(c)/Z^(0)), are like- 
wise obtained from ratios of partition functions Zq{c) with boundary conditions denoted by a <i 
dimensional vector c such that the q spin states = 0, I, ...q — I are cyclically shifted across the 
boundary, i.e. , = + c,- mod q, with = 0, 1 , . . . <7 — 1 . For q = 2, the Ising model, this 

simply means anti -periodic in direction /, if c, = 1 , and periodic otherwise. 

Interfaces(vortices) are suppressed in the low(high) temperature ordered(deconfinement) phase 
below(above) Tc. Here, the ratios /?g(c) = Zq{c)/Zq{0) and Rk(k) = Zj^(k) /Zi^{0) both tend to zero 
as we approach the thermodynamic limit. Complementary to that, in the disordered(confined) 
phase, it is the interface(vortex) free energies that tend to zero such that the ratios Rq, approach 
1 for all boundary conditions. Only at the critical temperature T = Tc do these ratios converge to 
non-trivial and universal fixed points. In 2 dimensions these universal numbers <Rc < 1 have 
been obtained exactly [9], in terms of Jacobi theta functions, for all cyclic boundary conditions 
c = (m, n) and all Potts models with 2"'' order transitions, i.e., for q = 2,3 and 4. 

These universal ratios can be used to determine the critical couplings jS^ of the transition in the 
2-1-1 dimensional gauge theories with high precision [10], by requiring that j3 = jSc at Rt = Re- For 



SU{2), where. 



R. 



1,0) 



1/(2^/4 + 1), and = (2^/4 -l)/(23/Vi: 



(2.1) 



from the 2d square Ising model, for example, this allowed us to determine the critical couplings 
for lattices with up to Nt =9 links in the Euclidean time direction with a typically by two orders 
of magnitude increased precision where there were literature values available [11, 12, 5]. The high 
precision allowed us to reliably determine the subleading 1 /A/^^-corrections to the linearly increasing 
behavior of j3c with A', near the continuum limit [10], 

pciNt) = 1.5028(21)A^,+0.705(21) -0.718(49) — . (2.2) 

Nt 

The slope of the leading term determines the critical temperature Tc in units of the dimensionful 
continuum coupling gj^ of the (2-i-l)<i theory, yielding Tc/g^ = 0.3757(5). By standard arguments, 
from this result one can also read off the temperature dependence of the coupling at a given A',, 

0.270(2) gj 



si 



Nr 



Tc l+t 



1 . 



(2.3) 



For each fixed Nt, we thus have precise control of the temperature by varying the lattice coupling 
j8. The physical length of the spatial volume then follows from L = aN^ = T^^{Ns/Nt) /{l+t). 
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With these we can perform a finite-size-scaling 
analysis near the 2"'^ order phase transition where gen- 
eralized couplings such as the vortex-ensemble ratios 
Ri^, for sufficiently large L, only depend on L^^^t, with 
the correlation-length critical exponent v, in a univer- 
sal way. In particular, for {2+l)d SU{2), this depen- 
dence is determined by that of the corresponding ra- 
tios Rq=2 of finite-volume partition functions of the 
2d square Ising model, with V = 1, which have been 
obtained exactly for all combinations of periodic and 
anti -periodic boundary conditions in [13]. 




Figure 2: Universal scaling function Fj{x) 
from the exact finite-volume partition func- 
tions of the N xN square Ising model in [13]. 



For one antiperiodic direction, the corresponding exact interface free energies per temperature 



InR. 



ap: 



where Rap stands for Rq=2{c] 



with either c = (1,0) or (0, 1) which are the same 
on an A'^ X A'^ square lattice, are plotted over the finite-size-scahng variable x = Nt for N = 100 
and = 1000 in Fig. 2. The asymptotic slope in the ordered phase {t < 0) for large A'^ approaches 
21n(H- \/2) which follows from Onsager's famous result [14] for the tension of a straight interface 
in the infinite volume limit, a/ = 2^ + lntanhiC' with spin-coupling per temperature K = J/T, near 
criticality at Kc = ln(l + \/2) /2. Deviations from that, e.g., observed for A^ = 100 around x = —5, 
are due to subleading (non-universal) finite-volume effects. 

We have calculated the SU{2) vortex-ensemble ratios Rk{k) with the algorithm of [15] as 
previously for the {'i+\)d theory in [16], here in {2+\)d for Nt = 2 to 10, each with spatial sizes up 
to Ns = 96 and in a suitable window of jS's around jSf to test the finite-size scaling (FSS). For each 
fixed A^, we generally observe very good scaling of the available data for all A'^, when plotted over 
the FSS waiiable X = LTct = {Ns/Nt)t{l5)/{l + t{l5)) (assuming v = 1 for the 2d Ising universality 
class) with f(j3) from Eq. (2.3). For the resulting vortex free energies (for k = (1,0) and (1,1)) 
we then obtain accurate one-parameter fits via Fk{x) = Fj{—Xx) to the exact universal scaling 
functions Fj{x) (for (1,0) and (1, 1) b.c.'s). This determines the single non-universal parameter A 
which is dimensionless and relates the SU (2) FSS variable to the Ising one (whereby the minus 
sign reflects the interchange of high and low temperature phases between the two). 



3. Self-Duality in 2+1 Dimensions 

The 2d ^-state Potts models are self-dual for all q. This has long been known for infinite 
systems. A particularly simple proof based on the random bond-cluster representation was given 
in [1]. Far less is known in a finite volume with translationally invariant boundary conditions, 
however, where ensembles with different boundary conditions mix under duality transformations. 
One case where this is known is the 3d Ising/Z2 -gauge theory system [17]. For the Potts model, 
using the random cluster methods developed in [9], we were able to obtain the following exact 
duality transformations for all g on a finite 2d torus as discrete 2d Fourier transforms over all 
ensembles Z^'"'"^ with cyclically shifted boundary conditions c = {m,n), 

Z^-)(^)= (^^)''^"^-Ief('™+-)zf")(^) m,n,r,s = 0,l,...q-l, (3.1) 
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Figure 3: The ratio of the partition function with twist k= (1,0) over the periodic ensemble compared to 
that of one unit of electric flux relative to the no-flux ensemble (left), Re ~ Zi,{e)/Ze{Q) with e = (1,0), and 
to its mirror image (right). The data here is obtained for Nt =4 and various spatial volumes up to Ns = 96. 

where A'^sites is the total number of sites of the graph on the 2d torus {i.e., on an N x N square 
lattice), K = J/T the coupling per temperature, and K its dual obtained from (e^ — 1 ) (e^ — l)=q, 
with temperature mirrored around criticality a.t K = K = Kc = ln{l -\- ^/q). 

For q = 2, and with K ^ K/2, conventionally, this relation reduces to an analogous known 
result for the Ising model [18]. For ^ = 3 it agrees with a corresponding result [19] obtained for the 
planar or vector Potts model which is equivalent to the standard one in that case. The proof of the 
general formula in Eq. (3. 1), for all <7-state Potts models on a 2d torus, will be given elsewhere [6]. 

Here we only note that the pattern in the finite-volume duality transforms, by the discrete 
Fourier transform over the ensembles with cyclically shifted boundary conditions, is precisely the 
same as that between the twisted vortex ensembles Zk(k,irt) in the d + I dimensional pure SU (N) 
gauge theory and those with fixed units of electric (and magnetic) flux Zg{e,m) [7], 

1 



Ze{e,m) 



Zk{k,m) 



(3.2) 



kez', 



obtained from a J-dimensional Zft^-Fourier transform over all temporally twisted ^-b.c.'s. We are 
not interested in magnetic flux and drop the argument in again, in the following. The role of the 
electric flux ensembles is best understood in terms of the translationally invariant flux between a 
fundamental color charge at some point x in the finite volume and its mirror charge at .x + eL in a 
neighboring volume in the direction of the flux ? [16], it contains no ultraviolet divergent perimeter 
and no short-distance Coulomb contributions [20] , 

Z,(?)/Z,(0) = ^{tr{Pix)P'{x + eL)))^^_^^^, (3.3) 

where the P's are untraced fundamental Polyakov loops, including any non-trivial transition func- 
tions in the time direction, and the subscript indicates that their expectation value is taken in the 
enlarged no-flux ensemble 2^(0) = Y^l^kik) /N"^ . Without loss, the form in (3.3) assumes periodic 
spatial b.c.'s and is manifestly invariant under spatially periodic gauge transformations. 

Of course, the duality relations provide exact maps between the high and low temperature 
phases in the spin models, valid at all temperatures and beyond universality. We do not have 
such exact duality relations for the gauge theories, but near criticality, where K ■(r^ K amounts to 
X -H- —X, we expect that the spin model dualities are reflected in analogously dual relations between 
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Figure 4: Vortex free energy Fk{x) in SUil) above r^. (left), and A(A^,) from fits to Fif{x) = Fi{-~Xx) (right). 



the behavior of the spatial center vortices and the confining electric fluxes. This is indeed the case 
in a rather wide universal scaling window around Tc and, in particular, for 2+\d gauge theories the 
self-duality of the spin models implies that the electric-flux free energies (e) are mirror images 
of the center vortex free energies, i.e., Fe{x) = Fj^{—x), as demonstrated for SU (2) in Fig. 3. 

4. Universal Amplitude Ratios and Continuum Limit 

The interface tension a/ in the ordered phase of the spin model corresponds to the dual string 
tension d for spatial center- vortex sheets which provides the dual area law for the spatial 't Hooft 
loops above Tc. Both vanish near the 2"'^ order phase transition, e.g., d ^ for t ^ 0+, with a 
critical exponent that is related to the coiTclation length critical exponent v via the hyperscaling 
relations a.?, ji = {d — \)v , i.e. jJ. = v ind = 2 spatial dimensions. The product of interface tension 
and correlation length (to the power J — 1) is an example of a universal amplitude ratio [21]. If 
we use the exponential correlation length in the disordered phase of the spin model, (§+p, then the 
universal constant given by Oi{^^^py^^ = Rag^p turns out to be exactly 1 for g = 2, 3 and probably 
also for <7 = 4 in 2 dimensions. In the gauge theory the exponential correlation length of Polyakov 
loops in the confined phase is given by string tension and temperature as (a/r)^'. Therefore, we 
have a universal relation between the string tension below and the dual string tension above Tc, 

^ = Ri,.,(^lTY-\ or a = pTc\-tr + -- and a = R+^^^ipTcY'-h^'-'^' + ■.■ , (4.1) 

where at leading order in the reduced temperature t the only unknown is the single unique constant 
p that we have introduced here. This constant can be extracted from the vortex free energies in the 
high temperature phase at asymptotically large x = TcLt^, where Ft{x) = R^^^^p'^'^^x^^^ ,hy direct 
numerical simulations. This was done in [16]. It is numerically very hard and expensive because 
larger and larger lattices are needed to stay within the universal scaling window as x is increased. 

This is demonstrated for {2+l)d SU (2) in Fig. 4 where the universal scaling function is known 
exactly from the Ising model. While its linear large-.x behavior, Fk{x) — )• px, is approached with 
increasing lattice size, even for the A^^ = 96 lattice it is still difficult to exti'act the asymptotic slope 
reliably from the data. But we know the asymptotic slope of the universal scaling function Fi{x) 
which is 21n(l + V2). Therefore, from Fjt(x) = Fj{—Xx) we obtain p = 2Aln(l + ^/2) with A 
from our one-parameter fits at small x where we have very accurate data. Moreover, we have 
determined this single non-universal parameter for A'^ = 4, 5 . . . 10 which allows a polynomial fit, 
^{Nt) = Xoo + b/Nf + c/Nf, with the extrapolated result Ac» = 1.354(25). This then determines the 
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leading behavior of the continuum string tension and its dual around the phase transition, 

a = Ar/21n(l + \/2)|f| + --- , f^O", and d = XT,.2ln{l + Vl) t + ■ ■ ■ , f^0+. (4.2) 

5. Conclusions 

Here we presented results for center vortex free energies in the 2+1 dimensional SU (2) gauge 
theory whose universal properties near the deconfinement transition can be analyzed in terms of 
many exact results from the 2d Ising model. We used these to precisely locate the transition, to 
demonstrate finite-size-scaling, and to extract the continuum string tension and its dual from univer- 
sal amplitude ratios together with one-parameter fits to exactly known universal scaling functions. 

We found an exact duality transformation for the ^-state Potts models on a finite 2d torus and 
demonstrated that their self-duality is reflected in the {2+l)d gauge theory: The free energies of 
spatial center vortices are mirror images around Tc of those of the confining electric fluxes. 

More results for SU (3) and SU (4) will be presented elsewhere, see also [6]. The relevance of 
center symmetry and center vortices for confinement in full QCD when including the electromag- 
netic interactions of the fractionally charged quarks are discussed in [22]. 
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